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THE  SOLUTION  OP  BOUNDARY  LAYER  EQUATIONS  BY  THE 
DIFFERENCE  METHOD 

I.  Yu.  Brailovskaya  and  L.  A.  Chydov 


In  this  article  we  present  the  methods  and  results  of  the  numeri¬ 
cal  solution  of  the  system  of  equations  of  a  two-dimensional  steady- 
state  boundary  layer  on  a  plate.  The  examined  method  can  be  applied 
in  the  particular  case  of  a  homogeneous  incompressible  fluid  and  in 
the  general  case  of  a  multi-component  compressible  gas  taking  into 
account  heat  transfer  and  diffusion. 

The  system  of  equations  of  the  boundary  layer  In  an  incompressible 
fluid,  as  we  know  [l],  has  the  form: 


..  i  „  1  do  ,  „  d*u 

u  — - — h  o - =  —  - --[-y - 

dx  dy  f  dt  dy* 

+  JL  „  o. 

dx  dy 


(1) 

(2) 


The  boundary  conditions,  taking  into  account  the  possible  passage 
of  the  fluid  across  the  surface  of  the  plate  (so-called  injection), 
are  set  as  follows: 


«  =  0,  o  =■  ob(at)  y  =»  o, 

u  =*  U(x)  *hm  y  =  00 
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(here  vw  (x)  is  the  transverse  velocity  of  the  fluid  on  the  plate 
surface,  differing  from  zero  when  injection  is  involved). 

In  the  case  of  a  compressible  gas  with  heat  transfer  and  diffusion, 
for  a  binary  mixture,  the  boundary  layer  equations  in  dimensionless 
form  (see  [2])  are 


On  tin  0  /  da  \  Op 

-ST1"1 

yt  1L  -vAe. JL  f-n  J_ 
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*  \  Ox  Oy  )  Oy  \  Oy  )  Pr.  r 
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(5) 

(6) 

(7) 

(8) 


The  boundary  conditions 


«  -  0;  o  =  tv.;  -  (I  -  c)  -Jf- ;  9  «  1 

n  —  1 ;  C  -----  0;  0  --  0  »htn  y  -■  co 


y  =  o. 


(9) 

(10) 


Here  we  use  the  designations  adopted  by  Avduyevskiy  and  Obroskova 
[2],  The  unknown  functions  are  u,  v,  c  and  0;  the  remaining  values 
are  known  functions  of  v,  c,  and  0,  or  are  constants. 

We  will  examine  Systems  (l)-(2)  and  (5)- (8)  with  corresponding 
boundary  conditions  (3),  (4),  and  (9)#  (10),  assuming  that  for  a 
certain  x  =  xQ  all  of  the  unknown  functions  are  given. 

For  the  numerical  solution  of  the  boundary  layer  equations  in  an 
incompressible  fluid,  the  method  of  Integral  equations  and  different 
variants  of  the  Galerkin  method  are  very  effective.  The  success  of 
these  methods  was  due  to  the  comparative  simplicity  of  System  (1),  (2) 
and  the  simple  construction  of  the  functions  describing  the  boundary 
layer  in  an  incompressible  fluid.  Attempts  to  use  similar  methods  to 
solve  systems  of  type  (4)-(8).#as  a  rule,  lead  to  unwieldy  non-standard 
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calculations  which  are  inconvenient  to  carry  out  on  high-speed  compu¬ 
ters. 

In  recenc  years  several  reports  have  been  published  concerning 
the  use  of  difference  methods  for  numerical  calculation  of  the  solution 
of  boundary  layer  equations  (see,  e.g.,  [>],  [4],  [8J,  [9J).  The 
application  of  difference  schemes  in  such  problems  :.a  complicated  by 
the  fact  that  the  boundary  layer  equations  are  peculiar  when  y  =  0, 
since  when  y  *  0  the  coefficients  of  the  derivatives  of  x  become  zero. 

As  a  result,  simple  evident  schemes  cannot  be  applied  close  to 
boundary  y  =  0  since  here  the  conditions  at  Ax,  associated  with  the 
requirement  that  the  difference  scheme  be  stable,  are  found  to  be 
unsatisfiable  in  practice.  In  an  earlier  article  [J>\  boundary  layer 
equations  of  the  Mlses  form  are  used.  The  above-indicated  peculiarity 
disappears,  but  another  difficulty  arises  since  the  second  derivative 
of  the  unknown  function  reverts  to  infinity  when  y  *  0.  This  prevents 
the  use  of  difference  schemes  with  small  values  of  jr.  Because  of  this, 
Mitchell  [j]  and  Baxter  [4]  used  special  methods  to  calculate  unknown 
functions  In  groups  lying  close  to  the  straight  line  y  0. 

In  view  of  the  unwieldiness  of  these  methods,  it  is  not  convenient 
to  use  them  in  more  general  cases  when  the  number  of  equations  Increases, 
for  example  In  the  case  of  a  multi-component  gas.  In  this  article  an 
Implicit  difference  scheme  Is  used  which  permits  calculation  of  the 
entire  boundary  layer  in  a  simple  uniform  way,  convenient  for  digital 
computers.  Implicit  difference  schemes  for  boundary  layer  calculation 
are  also  used  In  other  articles  [8]  and  [9]. 

The  system  of  boundary  layer  equations  In  any  case  consists  of  one 
first-order  equation  (continuity  equation)  and  several  second-order 
equations.  Each  of  these  equations  is  quasi-linear  (linear  with  respect 

-> 
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to  derivative)  relative  to  its  "own"  corresponding  function,  i.e., 
that  function  which  enters  the  equation  as  the  first  derivative  with 
respect  to  x  and  the  second  derivative  with  respect  to  The  equations 
are  connected  only  through  coefficients  which  depend  on  unknown  func¬ 
tions  and  also  on  their  first  derivatives  with  respect  to  jr  (as,  for 
example,  in  Eq.  (8)). 

Each  of  these  second-order  equations  can  be  expressed  in  the 
following  general  form: 


a 


f +t-f 

Ox  dy 


(11) 


The  boundary  conditions  can  also  be  expressed  in  general  form: 

df 


uoy-- !-P/-‘Y  y-o. 


/  -c--  const 


tj  ro. 


(12) 

(13) 


We  will  describe  the  difference  scheme  used  by  U3.  We  will 
Introduce  on  plane  (x,  y)  a  basic  rectangular  network  comprised  of 
points  with  coordinates 

X  =*  -V#  -I-  iii\x\  -J  rnSy  (m,  /i  =,  0, )  ,2, . . .). 

At  these  points,  which  we  will  call  integral  for  conciseness,  we  will 
calculate  the  values  of  the  functions  of  f.  satisfying  equations  of 
type  (11).  To  reduce  the  notation  we  will  introduce  the  designation 

f">  :  / (vj  •?-  «A.v,mAy).  (1*0 


In  addition  to  the  basic  network  we  also  need  an  auxiliary 
network  of  points  with  coordinates 

m  x  r»  V  u  =  ("I  -l-  */=)  A y 

«-^+(.+vj4x. 

At  these  "half-interger"  points  we  will  calculate  the  coefficients 
&,  b,  d,  and  £.  We  will  use  designations  analogous  to  (14).  A  little 


later  we  will  describe  a  specific  method  of  calculating  the  coefficients 
with  respect  to  the  value  of  the  functions  of  £  given  in  the  ©roups  of  the 
basic  network^  for  simplicity  of  explanation  we  will  thus  assume  that 
the  coefficients  a,  b,  d,  and  £  are  known. 

Differential  Equation  (11)  is  approximated  with  the  aid  of  a  two- 
layer  implicit  six-point  scheme: 


,!*•.  - ct\) -!•  (I - *) </•  „ 

.  *  %/m  .  —  ■  —  — — - . 

-1*  2  \g 


HmVl  ~  IV ')  ~  dm+-',  (  (/«  '  - 


(15) 


J.  n  _ A  l+l  //«)  rfw-,:,(/3i  —  /m-l)  , 

(%)«  • * 


where  ^  <  s  <  1. 

The  boundary  conditions  are  approximated  by  the  difference 
equations 

=  Y.  (l6) 

/.V,V7l,-c.  (17) 

Here  M(n  +  i)  is  the  number  of  the  last  group  point  for  which  unknown 
functions  in  the  layer  with  the  number  n  +  1  are  calculated.  Details 
on  the  selection  of  M(n  +  1)  will  be  described  below. 

The  calculation  will  be  conducted  layerwise,  consecutively  convert¬ 
ing  from  the  layer  with  the  number  n  to  the  layer  with  number  n  +  1. 

Let  us  assume  all  of  the  unknown  values  on  the  layer  with  the  number 
n  are  known.  Then,  if  the  values  of  the  coefficients  a,  b,  d,  and 
introduced  in  Difference  Equations  (15)  are  known.  Systems  (15),  (16), 
(17)  determine  the  values  of  £  in  the  layer  with  the  number  n  +  1. 

This  system  will  be  determined  by  the  usual  step  method  [5),  [6]. 
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We  will  find  step  coefficients  In  the  recursion  ratio 

fn  —  An  fn  »  Iin 
I  hi  --  -  'm/wiI  ;*  t 


by  known  formulas  consecutively  from  m  to  m  +  1.  Coefficients  AQ  and 

Bq  are  determined  with  the  aid  of  Boundary  Conditions  (16).  The  values 

of  fn  of  the  unknown  function  are  calculated  by  the  inverse  step  with 
m 

respect  to  Relation  (18)  with  the  aid  of  Boundary  Condition  (17). 

The  number  of  points  M(n  +  l)  in  the  layer  with  number  (n  +  l) 


is  determined  from  the  condition  that  the  difference  Bnouia 

have  a  sufficiently  small  absolute  value.  This  condition  provides  a 
smooth  transfer  of  £  to  its  limiting  value. 

In  practice,  M(n  +  l)  is  selected  as  follows.  First  we  will 
assume  M(n  +  l)  =  M*  =  M(n),  we  will  find  (as  has  been  described 
previously)  and  and  we  will  verify  the  condition  of  smooth 

conjugation: 

I  /.V  —  /"i'-il  < «.  (19) 

where  e  is  a  certain  small  positive  number.  If  (19)  is  not  fulfilled, 
we  will  assume  M  =  M"  =  M(n)  +  1  and  will  find  fjjn^  and  Thus, 

we  will  replace  the  missing  values  of  all  unknown  functions  in  the 
layer  with  number  n  with  their  limiting  values.  If  condition  (19) 
again  is  not  fulfilled,  we  will  Increase  M  to  unity,  and  so  forth, 
until  (19)  is  satisfied. 

After  we  obtain,  in  the  layer  with  number  (n  +  l),  the  values 
for  all  functions  of  £  satisfying  equations  of  type  (11 )  the  transverse 
velocity  v  is  found  from  the  difference  equation 

ihr  &■“)'»'  -  but]  -I-  - 1- !(,«)- :  a  -- 1  + 

(20) 

which  approximates  the  continuity  equation  (the  density  p  entering  into 
this  formula  depends  only  on  the  temperature  and  concentration  of  the 


should 
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components,  l.e. ,  on  functions  of  £,  and  therefore  we  now  can  consider 
p  as  a  known  value,  as  are  the  coefficients  a,  t>,  d,  and  £  of  the 
equations  of  type  (li) ) .  Efcr  Eq.  (20)  we  can  express  v^^2  in  terms 
of  v£+1/^2.  The  initial  value  of  v^+1/2  is  known  because  of  one  of 
Boundary  Conditions  (j).  Thus,  we  will  find  _v  in  the  layer  with  half¬ 
integer  number  n  +  1/2  at  the  points  with  whole  number  m. 

We  will  now  explain  the  method  of  calculating  the  coefficients 
entering  into  Difference  Equations  (15).  Obviously  from  Eqs.  (5)— (8) 
the  coefficients  a  and  £  depend  only  on  the  functions  of  £;  the  coef¬ 
ficients  £  and  £  can  depend  on  derivatives  of  the  type  df/dy.  For  the 
calculation  of  the  values 


a„  ‘  £>" +l* 

we  must  know  the  values  of  the  functions  of  £  and  p  in  the  groups  of 
the  auxiliary  network.  Using  linear  interpolation,  we  will  assume: 


/-,‘=Y(/r*+/fA 

I"'.  =  y  (/T  +  /DA'.). 


(21) 

(22) 


For  the  calculation  of  e^+1/^2  we  will  assume: 

m 


w,*.«  j-  JkjrA  u  •  JL  «*'.  -C-V  (23) 

\0V)M  2  A,,  Sg  •  '  ** 

Here,  naturally,  a  difficulty  arises;  for  calculation  of  the  functions 
of  £  in  the  layer  (n+l)  we  must  know  the  coefficients  which  themselves 
depend  on  unknown  values  of  the  functions  of  £.  This  difficulty  is 
surmounted  in  the  usual  manner  —  with  the  aid  of  so-called  enumerations 
or  iterations. 

In  the  first  approximation,  in  the  right-hand  part  of  (21),  (22) 
dnd  (25)  the  superscript  n  +  1  is  replaced  by  n  and  we  find  as  previ¬ 
ously  described,  the  first  approximation  for  £  in  the  layer  (n  +  l), 
and  v  in  the  layer  (n  +  i/2).  We  will  substitute  the  first  approxima- 


tion  for  £  in  place  of  the  values  of  £  in  the  layer  (n  +  l)  in  Formulas 
(21),  (22),  and  (23),  we  will  find  the  second  approximation  for  £  in 
layer  (n  +  l),  etc.  For  the  calculation  of  the  coefficients  a  and  b, 
we  must  know  v  in  addition  to  £  in  layer  (n  +  l/2).  For  the  calcula¬ 
tion  of  the  first  approximation  we  will  use  the  values  of  v  in  the 
next  known  layer,  i.e.,  in  the  layer  with  the  number  n  -  i/2;  in 
subsequent  approximations  we  can  use  the  previous  approximations  for  v. 

We  compiled  3-6  enumerations.  The  function  of  v  changed  most 
noticeably  in  the  enumerations. 

In  practical  calculations  we  used  a  scheme  which  varied  somewhat 
from  the  described.  The  basic  difference  was  the  application  of  a 
non-uniform  network. 

All  of  the  functions  which  describe  the  boundary  layer  undergo 
(with  fixed  x)  the  greatest  changes  with  small  values  y*  with  an 
increase  of  all  derivatives  decrease.  Thus  the  step  of  Ay  also 
increases  with  an  increase  of  jr.  In  each  layer  there  were  m^  points 
with  step  Ay^,  m^  points  with  step  Ay2  =  2Ay^,  and  m^  points  with 
step  Ay-j  =  2Ay2.  At  "junctions,"  i.e.,  at  the  points  m  »  m^  and 
m  =  mg,  first-  and  second-order  derivatives  were  approximated  with 
the  aid  of  the  values  of  the  functions  at  the  points  with  numbers 
m  -  2,  m,  m  +  1  (second-order  approximations  relative  to  Ay).  We  will 
not  present  here  the  obvious  changes  which  must  be  made  in  Difference 
Equations  (15 ),  (20)  and  recursion  formulas  for  the  step  coefficients 
A^,  Bffl  for  points  at  which  the  step  with  respect  to  Ay  changes.  The 
use  of  the  variable  step  with  respect  to  ^  made  it  possible  to  consid¬ 
erably  reduce  the  total  number  of  calculation  points.  At  "junctions" 
we  observed  no  disruption  of  the  smoothness  of  the  functions  to  be 
calculated.  The  width  of  the  boundary  layer,  generally  speaking. 
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increases  with  an  increase  of  x. 

In  the  course  of  the  calculations  we  changed  the  "averaging  param¬ 
eter"  appearing  in  Difference  Equation  (15),  taking  into  account 
the  nature  of  the  change  of  the  functions  of  £  with  respect  to  x. 

As  is  obvious  from  Eqs.  (l),  (2)  or  (5)-(8)»  the  transverse  veloc¬ 
ity  v,  when  x  =  xQ,  is  determined  by  giving  values  of  all  remaining 
functions  when  x  *  xQ.  Thus,  for  example,  from  (l)  and  (2)  we  will 
have: 


Ov  .  du  1  Bp  ,  d*u 

a - \-  v - - - .j.  v-— 

Oy  Oy  f  d*  dy* 


(24) 


This  ratio  represents  a  usual  linear  differential  equation  relative 
to  v(xQ,  y).  Boundary  Conditions  (3)  determine  the  initial  value  of 
v(xQ,  0).  Numerical  Solution  of  Eq.  (24)  is  difficult  because  the 
coefficient  when  ^  reverts  to  zero  when  y  ■=  0. 

The  initial  function  v(xQ,  y)  can  be  obtained  with  the  aid  of 
the  iterative  process  shown  on  p.  7,  This  process  will  give  us 
v(xq  +  ,  y);  for  sufficiently  small  Ax  we  can  identify  this  function 

with  v(xQ,  y).  • 

When  calculating  the  initial  function  v(xq,  y)  we  noticed  a  great 
stability  of  the  Iterative  process  relative  to  perturbation  of  the  zero 
approximation  for  Convergence  of  the  iterative  process  was  not 
disrupted  even  for  very  large  deviations  of  the  zero  approximation  from 
the  true  function  v(xQ,  y).  The  error  decreased  with  the  rate  of  a 
geometric  progression  having  an  Index  close  to  1/2. 

The  stability  of  Difference  Equation  (15)  was  Investigated  in  a 
case  of  constant  coefficients  a,  ]3,  and  £  with  constant  steps  Ax  and 
Ay  of  the  network,  and  with  uniform  boundary  conditions.  The  invests 
gation  was  easily  conducted  by  means  of  the  Fourier  transform  [7J. 
Stability  occurs  with  an  arbitrary  convergence  of  Ax  and  Ay  to  zero. 
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Judging  from  the  results  of  practical  calculations,  the  facts  distin¬ 
guishing  the  actual  problem  from  this  model  case  (boundary  conditions, 
variable  step  with  respect  to  £,  presence  of  non-linear  relations 
between  equations  of  the  system)  do  not  disturb  the  stability  of  the 
scheme . 

The  error  of  the  method  presented  above  was  verified  by  comparison 
of  results  obtained  using  various  Steps  of  Ay  and  Ax,  and  by  comparison 
with  knowi  solutions  found  by  other  methods. 

For  a  comparison  with  the  familiar  solution  of  the  Blaslus  equa¬ 
tion  by  the  described  difference  method,  we  solved  System  of  Equations 
(l)-(4)  with  a  zero  pressure  gradient  and  zero  injection  velocity  vw. 
Calculations  were  made  from  Xq  =  O.Olj  the  solution  of  the  Blaslus 
equation  [1]  was  used  as  the  initial  data.  The  initial  layer  contained 
16  nodal  points,  not  counting  the  points  on  the  wall,  wherein  m^  «  10, 
mg  =  4,  m^  =  2j  the  values  of  the  corresponding  steps  were  Ay^  «  0.04, 
Ay2  =  0.08,  Ay^  =  0.l6,  and  Ax  =  0.005. 

Five  steps  according  to  x  were  made:  the  obtained  results  of 
function  u  did  not  differ  from  the  solution  of  the  Blasisus  equation 
by  more  than  0.25#. 

System  of  Equations  (5)-(10),  where  =  0,  was  solved  by  the 
above-mentioned  difference  method  for  several  variants. 

As  the  initial  data  in  all  variants,  we  used  the  solution  obtained 
in  a  previous  article  [2]  for  a  particular  case  of  boundary  conditions 
that  permitted  a  self-similar  solution. 

The  calculations  were  carried  out  from  Xq  =  0.01. 

1.  For  the  segment  0.01  £  x  <  O.O725  calculation  was  conducted 
with  self-similar  boundary  conditions  [2]  with  a  step  of  Ax  =  0.0025 
according  to  x.  In  the  initial  layer  18  points  were  taken,  not  counting 
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the  points  on  the  wall,  wherein  »  10,  mg  -  5>  Bj  ■  5j  the  steps 
according  to  £  were  Ay^  *  0.05j  Ay2  ■  0.1,  Ay^  -  0.2. 

With  an  Increase  of  x*  the  width  of  the  layer  Increased  In  the 
manner  described  above;  the  conditions  of  a  smooth  conjugation  on 
the  upper  boundary  of  the  layer  were  verified  only  for  the  concentra¬ 
tion  function  _c,  since  the  remaining  functions  go  to  the  asymptote 
when  y-»  »  mare  rapidly  than  £. 

The  functions  0,  and  u,  being  a  solution  of  System  (5)-(10) 
in  the  self-similar  case,  depend  only  on  a  certain  variable  similitude 
€  *  £(x,  y)*  Consequently,  the  values  0,  _c  and  u  obtained  by  the 
difference  method  with  various  values  of  x  (in  various  layers  with 
respect  to  x)  3hould  agree  during  conversion  to  this  variable,  if  our 
solution  agrees  with  the  exact  solution  of  System  of  Differential 
Equations  (5)-(10). 

From  the  exactness  of  the  coincidence  we  can  judge  the  error 
of  the  method.  An  examination  of  the  calculation  results  shows  that 
on  the  segment  0.03  <  x  <  0.0725  "scattering"  of  the  values  of  func¬ 
tions  for  identical  ?  does  not  exceed  2.5#  for  u,  3#  for  £,  and  for 
0,  for  the  values  of  Q  in  the  interval  1  >  6  >  0.5  "scattering"  does 
not  exceed  3#.  (The  indicated  values  of  0  belong  to  the  interval 
0  <  y  <  0.46  where  x  =  0.01;  the  total  thickness  of  the  layer  x,  in 
this  case  i3  equal  to  1.6.) 

The  difference  from  the  solution  obtained  in  a  previous  article 
[2]  did  not  exceed  3.5#  for  u,  4$6  far  c_,  and  10#  for  0.  The  change 
of  the  averaging  parameter  did  not  noticeably  influence  the  calcula¬ 
tion  results. 

2.  With  the  same  distribution  of  points  in  the  initial  layer 
as  in  case  1,  a  calculation  was  conducted  with  Boundary  Conditions  (9)# 


(10)  where  we  assumed  vw  *  const  (these  boundary  conditions  are  not 
self-similar).  The  step  according  to  x  was  taken  as  constant  and 
equal  to  0.0025.  The  calculation  was  conducted  up  to  x  «  0.1625.  The 
width  of  the  boundary  layer  was  increased  from  L*  1.6,  when  x  ■  Xq  « 

=  0.01,  to  L  =  7.2,  when  x  =  O.165.  The  value  of  L  in  this  and  In  all 
other  variants  was  determined  from  the  condition  of  smooth  conjugation 
for  the  concentration  function  £  of  the  upper  boundary  of  the  boundary 
layer. 

The  first  six  layers  according  to  x  were  calculated  with  six 
Iterations  and  s  =  1  in  order  to  acquire  rapid  attenuation  of  the 
initial  error.  Then  we  assumed  that  x  =  0.07  and  the  calculations 
were  carried  out  with  three  iterations.  The  calculation  results  are 
presented  in  the  table. 

3.  The  third  variant  differs  from  the  second  in  that  with  a 
doubling  of  the  boundary  layer  thickness  in  comparison  with  the  original 
the  steps  according  to  £  doubled.  At  the  same  time  the  step  according 
to  x.  doubled.  This  procedure  was  repeated  for  each  recurrent  doubling 
of  the  thickness  of  a  boundary  layer.  The  initial  step  according  to 
x  equaled  0.0025.  In  this  variant  the  calculations  were  carried  to 
x  =  0.245. 

A  comparison  of  variants  2  and  3  showed  that  the  difference  in 
the  values  of  u  did  not  exceed  5#  for  all  x;  the  same  was  true  for  v. 

The  difference  in  the  values  of  £  for  all  x  does  not  exceed  7%  for 
values  of  this  function  larger  than  0.1.  The  difference  in  the  values 
of  0  for  all  x  does  not  exceed  7%  for  values  of  this  function  larger 
than  0.3. 

For  small  values  of  £  and  9  pertaining  to  large  values  of  £,  the  dif¬ 
ference  in  the  results  of  these  variants  increases.  The  reason  obviously  is 
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that  the  step  according  to  £  far  from  the  wall  is  too  great  in  variant  J 

In  a  direction  from  the  upper  boundary  of  the  layer  to  the  wall 
the  error  rapidly  diminishes  as  is  obvious  f"om  a  comparison  of  the 
results  of  these  variants.  We  can  Judge  the  rate  of  diminishing 
according  to  the  value  of  the  step  coefficients  since,  as  is 
obvious  from  (l8),  the  introduction  of  error  e  into  the  values  £  and 
0  when  m  *  raQ  causes  the  appearance  of  an  error  of  the  order 

"b 

n  A,,e  in  these  values  when  m  -  nu  -  r. 
k-m0-5* 

In  our  calculations,  in  the  examined  region  far  from  the  wall  k^ 
changed  from  0.016  to  0.08  for  the  function  0  and  from  0.012  to  0.06 
for  the  function  £. 

For  estimations  of  the  influence  of  the  averaging  parameter  £ 
in  the  third  variant  after  the  sixth  layer  we  assumed  js  equal  to  l/2. 

The  results  obtained  in  this  case  (the  calculation  was  carried  to 
x  -  0.18)  differed  little  from  the  corresponding  results  of  the  third 
variant. 

4.  in  order  to  verify  the  accuracy  of  choice  of  step  Ay  in  the 
initial  phase  of  the  calculation,  control  calculation  with  a  more 
dense  network  (according  to  £)  was  carried  out  in  the  segment  0.01  < 

£  X  £  0.02. 

In  the  initial  layer  56  points  were  taken,  with  m^  *  20;  mg  =  10; 
m-j  «  6;  the  steps  according  to  £  were  taken  equal  to  0.025,  0.05,  and 
0.1,  respectively.  The  step  according  to  x  was  equal  to  0.0025.  The 
calculation  was  carried  out  rith  six  iterations  where  s  -  1.  The 
results  obtained  in  this  case  for  the  functions  0,  £,  u,  and  £  differ 
very  insignificantly  from  the  results  of  variants  2  and  >  at  the 
corresponding  points  of  the  network. 


Below  is  the  table  of  the  calculation  results  for  variant  2. 

In  the  calculation  we  used  the  following  formulas  for  determining 


the  coefficients  of  the  equations,  depending  on  u,  £, 

X  =  {-2,01c 1 4- 0.0203  _ 1 .2] . .  Q  r 

=  I0.9c>  (•  8,0c’  4  0,912c  j-  0,031  "l. 2  -0.5  9  '' 

p  38  (i  -O.SOHl  J-  I3,5<j* 

""  ,+3SSr^ ) 

c,  =  (i--O.50)#>»(l3ki5c  +  l), 

Drt  =  (l-0,50)‘.» 

•  cP.  -  cp,  =  1 3,15  (1  —  0.5  G)°.'» 

ej>.» 

k  =  0.SIG326,  Prfivc  =  0,151,  PrK=0.73. 


and  8: 

1 .5 


The  derived  formulas  pertain  to  the  case  of  a  hydrogen-air  mixture 


I 


*=6,012$ 


9 

« 

c 

II 

0.0 

1,00000 

0,41329 

0.00000 

4,2880 

0.1 

0,04315 

0>;6$8 

0,21135 

4,0584 

0.2 

0,87034 

0,20323 

0,44293 

5,3902 

0,3 

0,77103 

0.13539 

0.66056 

8,1549 

0,4 

0,03143 

0,05954 

0,82542 

11,375 

0.5 

0/.S479 

0,05980 

0.92265 

14.123 

0.6 

0,29527 

0.0 1024 

0.96970 

15,903 

0.7 

0,16152 

0.02722 

0,93878 

17,093 

0.8 

0 ,07635 

0,01825 

0,99634 

18,052 

0,9 

0,03309 

0,01207 

0,99376 

19,095 

1.0 

0,01509 

0,00787 

0,99939 

20.131  ‘ 

1.2 

0,00306 

0,00312 

0.99989 

21 .816 

1,4 

0.00058 

0,00113 

0,99998 

22,854 

1.6 

0,00011 

0,00036 

0,99999 

23,345 

1.8 

0.00000 

0.00000 

1.00000 

23,445 

' 

*=0,0425 

0.0 

I.CGOOO 

0.62526 

0.00000 

4,2880 

0.1 

0,98164 

0,51688 

0,06373 

3,9923 

0.2 

0,95982 

0,-16682 

0,13925 

3,9570 

0,3 

0,93101 

0.39007 

0,22565 

4,2393 

0.4 

0.90365 

0,32048 

0,32077 

4,8654 

0.5 

0,86818 

0,26016 

0.42133 

5.8209 

0,6 

0.S2699 

0,20965 

0,52287 

7,0488 

0.7 

0.77914 

0,16837 

0,62039 

8,4491 

0.8 

0,72192 

0,13517 

0 .70092 

9,9079 

0,9 

0.66308 

0,10877 

0.78138 

11,314 

1.0 

0.59119 

0,08790 

0.84411 

12,572 

1.2 

0.-14257 

0,05764 

0,92138 

14.442 

1,4 

0.29363 

0,03806 

0,96614 

15,663 

1.6 

0,17181 

0.02515 

0,98582 

16,482 

1.8 

0.0S873 

0,01650 

0,99439 

17,116 

2.0 

0.01084 

0,01066 

0,99789 

17,666 

2.2 

0.01702 

0,00673 

0.99924 

18,147 

2.4 

0 ,00653 

0,00414 

0,99974 

18.545 

2.6 

0.00231 

0,00244 

0,99991 

18,851  • 

2.8 

0.00078 

0,00136 

0,99997 

19,072 

3.0 

0,00024 

0,00069 

0.99999 

19,222 

3.2 

0,00005 

0.00028 

0,99999 

19,319 

3.4 

0.00000 

0.00000 

1.00000 

19,339 

x=0,0725 

0.0 

1 ,00000 

0.76164 

0,00000 

4,2880 

0.1 

0,99101 

0,71081 

0,02948 

4,1001 

0.2 

0.98062 

0,65399 

0,06506 

3,9771 

0.3 

0.90301 

0,39292' 

0,10721 

3,9444 

0.4 

0,95178 

0.52978 

0,15606 

4,0273 
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**-0,0725 


V 

« 

c 

u 

0 

0,5 

0,93590 

0,46695 

0,21140 

4,2474 

0.6 

0.92074 

0,10667 

0,27262 

4,6231 

0.7 

0,90004 

0,35061 

0,33873 

5.1536 

0.8 

0,87656 

0.29C38 

0,4081.3 

5.S314 

0.9 

0.85004 

0,25498 

0,48018 

6,6.382 

1.0 

0,82022 

0,21595 

0,55229 

7.5164 

1.2 

0.74896 

0,15392 

0,68717 

9,1795 

1.4 

0.G6055 

0,10954 

0,79820 

11,322 

1.6 

0,55544 

0,07325 

0,87907 

12,854 

1.8 

0,43933 

0,05620 

0,93216 

11.002 

2.0 

0.32295 

0,04055 

0,92!’  ' 

14,808 

2.2 

0,2)861 

0,02931 

0.2-SvJ 

15.369 

2.4 

0.13552 

0,02113 

0,93113 

15,788 

2.6 

0,07634 

0,01514 

0,99609 

16.137 

2.8 

0.03935 

0,01072 

0.S9S29 

16,453 

3.0 

0,01913 

0.00747 

0,99323 

16,7-11 

3.2 

0.00857 

0,00511 

0,939/1 

16,997 

3.4 

C. 00359 

0, 00340 

0,99959 

17.214 

3,6 

O.COl  12 

0,60219 

0.99936 

17,391 

3.8 

0,00052 

0,00135 

0.CS998 

17.528 

4.0 

0,00018 

0,60377 

0,99339 

17.529 

4.2 

0,00005 

0,00038 

0,00999 

17,697 

4.4 

0,00000 

0.00012 

0,90399 

17,757 

4.6 

0.00000 

O.CGOOO 

*=0,1025 

1,00000 

17.SQ6 

0.0 

1. 00000 

0,85698 

0.00000 

4,2880 

0,1 

0.99493 

0,82618 

0.01478 

4,1835 

0.2 

0.93914 

0,79015 

0,03276 

4.1004 

0.3 

0,98255 

0,74895 

0,05139 

4,0184 

0.4 

0,97505 

0.70302 

0,08012 

4,0389 

0.5 

0,96654 

0,65315 

0,11024 

4.0852 

0.6 

0,95683 

0,60019 

0.11501 

4.2028 

0,7 

0,94597 

0.51641 

0,18146 

4.4021 

0.8 

0,93367 

0,49231 

0,22317 

4,6930 

0.9 

0.91985 

0,43955 

0,27674 

5,0823 

1.0 

0,90137 

0,33929 

0.32873 

5.5719 

1.2 

0,86775 

0,29988 

0,44065 

6,8526 

1.4 

0,82242 

0,22711 

0,55673 

8,4198 

1.6 

0.70674 

0,17013 

0,66774 

10,1164 

1.8 

0.69913 

0,12751 

0,76193 

11,7634 

2.0 

0.61885 

0,09548 

0,84285 

.13,2149 

2.2 

0,52731 

0,07174 

0,90012 

14,3921 

2.4 

0.42S79 

0,05412 

0,93996 

15,2842 

2.6 

0,33010 

0.01C95 

0,96545 

15,9294 

2.8 

0.23S96 

0,03103 

0,95096 

16,3902 

3.0 

0,16181 

0,02349 

0,98994 

16.7313 

♦ 


I 


i 


16 


x^0tl02S 


* 


p 


* 


» 

• 

e 

II 

0 

3.2 

0.10223 

0,01772 

0,09189 

17,005 

3.4 

0.06021 

0,01327 

0,99750 

17,214 

3.6 

0.03311 

0,00985 

0,99882 

17,163 

3.8 

0,01706 

0,00721 

0,99916 

17,666 

4.0 

0,00828 

0,00519 

0,99976 

17,851 

4.2 

0,00380 

0,00367 

0,99990 

18,012 

4.4 

0,00166 

0,00253 

0,59996 

18.148 

4.6 

0,00068 

0.00168 

0,99998 

18.259 

4.8 

0,00027 

0,00106 

0,95999 

18,346 

5.0 

0,00009 

0.00062 

0,99999 

13,413 

5.2 

0,00003 

0,00028 

0,99999 

13,456 

5,4 

0.00000 

O.OoOOO 
x— 0,1325 

1.00000 

18,442 

.  0.0 

1,00000 

0,92620- 

0,00000 

4.2830 

0.1 

0.99716 

0.91022 

0,00082 

4,2:01 

0.2 

0,99301 

0 .89099 

0,01514 

4,2035 

0.3 

0,99032 

.  0,86315 

0,02524 

4.1321 

0.4 

0.9SG22 

0,81)13 

0.03713 

4,1807 

0.5 

0,93160 

0,81063 

0,05202 

4,205$ 

0.6 

0.97010 

0,77575 

0,06931 

4,2650 

0.1 

0.97055 

0.73693 

0,08961 

4.3661 

0.8 

0,96397 

0.69473 

0,11317 

4.5176 

0.9 

0.95660 

0,64964 

0,11017 

4,7287 

1.0 

0,91835 

0,60251 

0,17072 

5.0077 

1.2 

0.92888 

0,50643 

0,24238 

5.3158 

1.4 

0.90184 

0.41397 

0,32683 

6,9524 

1.6 

0.87511 

0,33093 

0,42093 

8.3978 

1.8 

0,83966 

0,20033 

0,51995 

10.062 

2.0 

0,79653 

0,20272 

0,61805 

11,893 

2.2 

0.711S0 

0.15702 

0,70912 

13,698 

2.4 

0,68317 

0,12144 

0,78817 

15,370 

2.6 

0,61223 

0.09401 

0,85238 

16,818 

2.8 

0.53217 

0,07296 

0,90141 

17.998 

3.0 

0,41615 

0,05679 

0,93675 

18,913 

3.2 

0,35861 

0,04432 

0,96095 

19,592 

3.4 

0.27199 

0,03165 

0.97677 

20.085 

3.6 

0,20023 

0,02708 

0,93666 

20,443 

3.8 

0.13793 

0,02113 

0,99259 

20,715 

4.0 

0,059-31 

0,01612 

0,99603 

20,935 

4.2 

0,05523 

0,01268 

0,99794  • 

21,126 

4.4 

0.03209 

0,00972 

0,99896 

21,301 

4.6 

0,01766 

0,00736 

0,99919 

21,461 

4.8 

0.00921 

0,00551 

0,03976 

21.608 

5.0 

0,00161 

0,00105 

0,59989 

21,739 

5.2 

0.00219 

0,00292 

0,99995 

21,854 

5.4 

0.C0100 

0,00205 

0,99997 

21,951 
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jr =0,1325 


y 

1 

C 

II 

0 

5.6 

0,00014 

0,00138 

0,99099 

22.031 

5,8 

0,00018 

0,00088 

0,99999 

22,095 

6.0 

0,00006 

0,00050 

0,09999 

22.145 

6,2 

0,00002 

0,00022 

0,99999 

22,185 

6,4 

0.00000 

0.00000 

*=0,1625 

1,00000 

22,187 

0.0 

1.00000 

0,97444 

O.OOCOO 

4,2880 

0.1 

0,90867 

0,96888 

0,00222 

4,2750 

0-.2 

0,99717 

0.96208 

0,t049l 

4,2691 

0.3 

0,99549 

0,95381 

0,00326 

4,2721 

0,4 

0,99361 

0,94379 

0,01231 

4,2859 

0,5 

0,99149 

0,93177 

0,01722 

4,3134 

0,6 

0,98912 

0,91744 

0,02314 

4,3580 

0.7 

0.98645 

0,90051 

0,03028 

4,4233 

0.8 

0,03346 

0,88072 

0.03832 

4.5139 

0,9 

0.9S012 

0,85782 

0.04S95 

4,5352 

1.0 

0,97637 

0.83166 

0,00094 

,4,75  >2 

1.2 

0.96751 

0.76942 

0,09136 

5,2  7,1 

1.4 

0,95618 

0,69511 

0,13176 

5.2634 

1.6 

0,91286 

0,61210 

0,18523 

G,„7IS 

1.8 

0.90616 

0,52554 

0,21615 

8.3209 

2.0 

0.90386 

0,11096 

0,31942 

10X52 

2.2 

0.8S141 

0,36293 

0,10090 

12,113 

2.4 

0,85213 

0,29432 

0,48735 

14.513 

2.6 

0.81729 

0,23619 

0,57475 

17,068 

2.8 

0,77604 

0.18829 

0.65876 

19,674 

3.0 

0,72751 

0,14961 

0,73544 

22,197 

3.2 

0.67108 

0,11875 

0,80188 

24,519 

3.4 

0,60675 

0,09434 

0,85664 

25.555 

3.6 

0,53542 

•  0,07507 

0,89967 

28.262 

3.8 

0,45915 

0,05987 

0,93202 

29,634 

4.0 

0,38103 

0,04784 

0,95537 

30,697 

4.2 

0.30481 

0,03827 

0,97158 

31 , 195 

4.4 

0,23425 

0,03053 

0,98244 

32,081 

4.6 

0.17248 

0,02448 

0.9S916 

32,509 

4.8 

0.12144 

0,01953 

0,99335 

32.827 

5.0 

0.08167 

0,01551 

0.99651 

33,073 

5.2 

0.05247 

0,01225 

0.99S07 

33.276 

5.4 

0,03223 

0,00960 

0,99596 

33,452 

5.6 

0.01895 

0,00746 

0,93916 

33,611 

5.8 

0.01068 

0,00573 

0,99973 

33.756 

6.0 

0,00579 

0,00434 

0.99986 

33,887 

6.2 

0,00302 

0,00324 

0,99993 

34,004 

6.4 

0,00152 

0,00236 

0,99997  • 

34,106 

6,6 
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